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ABSTRACT 

We report on 5 Chandra observations of the X-ray afterglow of the Gamma-Ray Burst GRB 060729 
performed between 2007 March and 2008 May. In all five observations the afterglow is clearly detected. 
The last Chandra pointing was performed on 2008-May-04, 642 days after the burst - the latest 
detection of a GRB X-ray afterglow ever. A reanalysis of the Swift XRT light curve together with 
the three detections by Chandra in 2007 reveals a break at ^1.0 Ms after the burst with a slight 
steepening of the decay slope from a = 1.32 to 1.61. This break coincides with a significant hardening 
of the X-ray spectrum, consistent with a cooling break in the wind medium scenario, in which the 
cooling frequency of the afterglow crosses the X-ray band. The last two Chandra observations in 2007 
December and 2008 May provide evidence for another break at about one year after the burst. If 
interpreted as a jet break, this late-time break implies a jet half opening angle of ^ 14° for a wind 
medium. Alternatively, this final break may have a spectral origin, in which case no jet break has been 
observed and the half-opening angle of the jet of GRB 060729 must be larger than ^ 15° for a wind 
medium. We compare the X-ray afterglow of GRB 060729 in a wind environment with other bright 
X-ray afterglows, in particular GRBs 061121 and 080319B, and discuss why the X-ray afterglow of 
GRB 060729 is such an exceptionally long- lasting event. 

Subject headings: gamma rays: burst. X-rays: burst. X-rays: individual (GRB 060729), Swift, Chandra 



1. INTRODUCTION 

Gamma-Ray Bursts (GRBs) are the most energetic 
transient events in the Universe. The most accepted 
mode l to explaiii this phenomenon is the "firebal l" model 
(e.g. iMeszaros 1 120061 : IZhmg fc Meszaros1[200l and the 
progenitor of long GRBs are believed to be massive stars 
on the order of 30 solar masses or more. The isotropic 
energies inferred from the observed fluences are often of 
order 10^^ — 10^^ ergs. On the other hand, typical super- 
nova explosions are of the order of only 10^^ ergs. One 
way to solve this "energy problem" is by assuming that 
the radiation is collimated into a jet. One prediction 
of the "fireball" model is that the afterglow decay rate 
increases when the relativistic beaming angle equals or 
exceeds the physical jet opening angle as the jet decel- 
erates in the su rrounding medium (e.g. iRhoadsl [19991 : 
ISari et al.l[l999h . This can be seen as an achromatic jet 
break in the light curve, with a typical decay slope after 
the jet break of order a = 2 (e.g. Zhang et al. 2006; 
Nousek et al. 2006). The time of this jet break can be 
used to infer the jet opening angle by e.g. the relation 
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given by iFrail et all (j2001| ) . The measurement of the jet 
break time is therefore most critical for understanding 
the energetics of GRBs. 

Before the la unch of the Swift Ga mma-Ray Burst Ex- 
plorer mission (jGehrels et al.|[20M) . putative jet breaks 
were found at optical or radio wav elengths, typicall y 
a few days after the burst (e.g. iFrail et aD 1200 ll ). 
Since its launch, Swift has detected roughly 400 bursts 
(end of 2008) and typically observes them for up to 
one or two weeks after the trigger. For the major- 
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reason could be that jet breaks occur in X-rays at 
much later times than previously thought and the follow- 
up observations by Swift are not late enough in time 
to de tect a jet break , as st udies by iWillingale et al.l 
(120071) and iSato et all (I2007D suggest. As an alterna- 
tive. iCurran et al.l ( 2008| ) suggested that the jet breaks 
are hidden and the light curves are mis-interpreted as 
a single power law decay although a jet break is there. 
Only a handful of afterglows have been followed for more 
than a month by Swift, because usually the X-ray af- 
terglow fades below the Swift X-Ray telescope (XRT; 
IBurrows et al.l [20051) detection limit ^ 10~^^ ergs s~^ 
cm~^ roughly a week after the trigger. One of these ex- 
ceptions is the bright X-ray afterglow of GRB 060729 
which was detected by Swif t XRT even 125 days after 
the burst (jCrupe et al.ll2007D . 

Th e Swift Burst Alert Telescope (BAT; iBarthelmvl 
[SMI triggered on GRB 060729 on 2006 July 29, 19:12:29 
UT (IGrupe et al.l [20071) and a re dshift of z = 0.54 
was measured (|Thoene et alll2006f). The XRT and th e 
UV/Optical Telescope fUVOT: iRoming et al.l [2005[) . 
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started observing the burst about 2 minutes after the 
trigger. The Swift UVOT was able to follow this after- 
glow in the UVWl filter up to 31 days after the BAT trig- 
ger. In X-rays Swift^s XRT was still detecting the X-ray 
afterglow of GRB 060729 at the end of 2006 November, 
125 days after the burst (jGrupe et al.ll200"7i ). However, 
by 2006 December the Swift-XRT detection limit was 
reached and only a 3a upper limit could be given for the 
63.5 ks exposure time obtained in December 2006. By 
that time the X-ray afterglow of GRB 060729 did not 
show any clear evidence for a jet break, gi ving a lower 
limit on the jet opening angle of 6* = 28° (IGrupe et al.l 
I2007D based on the assumption of a constant circumburst 
medium. In order to extend the light curve of this ex- 
ceptional X-ray afterglow, we observed it five times with 
Chandra ACIS in 2007 and 2008. 

We report on the detections o f the X-ray afterglow of 
GRB 060729 (|Grupe et al.|[2007[ ) with Chandra up to 642 
days after the burst - the latest detection ever of an X-ray 
afterglow of a GRB at cosmological distance. Previously 
the burst with th e latest detection of an X -ray afterglow 
was GRB 030329 (jTiengo et al.|[200l[200l . which had a 
detection 258 days after the burst by XMM-Newton. Our 
paper is organized as follows: in §[2]the observations and 
data reduction are explained, in §[3] the measurements of 
the X-ray light curve are shown, and in §S] we discuss 
the implications of this light curve. Throughout this pa- 
per the X-ray flux dependence on time and frequency is 
defined as F oc t~'^v~^ . Luminosities are calculated as- 
suming a ACDM cosmology with fiivi = 0.30, Ha = 0.70 
and a Hubble constant of = 71 km s^^ Mpc^^ cor- 
responding to a luminosity distance = 3064 Mpc for 
GRB 060729. 

2. OBSERVATIONS AND DATA REDUCTION 

Chandra observed GRB 060729 three times between 
2007 March 16 and 2007 June 30. Two very late-time 
Chandra observations were performed in 2007 Decem- 
ber/2008 January for 72.7 ks and in 2008 April/May for 
117.3 ks. Due to pitch angle constraints some of these 
had to be split into several visits. All observations, with 
start and end times and exposure times, are listed in 
TablelH 

All of these observations were performed with the stan- 
dard 3.2s readout time in Very Faint mode on the on- 
axis position on the back-illuminated ACIS-S3 CCD. 
Data reduction was performed with the Chandra analysis 
software CIAO version 4.0 and the calibration database 
CALDB version 3.4.3. In order to reduce the ACIS par- 
ticle background all Chandra stage 1 event data were 
reprocessed using CIAO acis_process_events with the 
VF mode cleaning. Only ACIS grades 0, 2, 3, 4, and 6 
were selected for further analysis. The background was 
further reduced by using only photons in the 0.5 - 8.0 keV 
energy range. Before further analysis, the observations 
were combined into one event file each using the CIAO 
task merge-all. Source photons were selected in a circle 
with a radius r = 1 and background photons in a close- 
by source-free region with a radius r = 10 . Count rates 
were converted into fluxes by using PIMMS version 3.9b 
using the parameters from the spectral fits to the Swift 
data after the break at 1 Ms after the burst (see below) 
with an absorption column density iVn = 1.34 x 10^^ 



cm^^ and an X-ray spectral index /3x= 0.89. 

A description of the reduction and analy sis of the Swift 
data can be found in iGrupe et al.l ()2007D . For display 
purposes and fitting the late-time light curve we rebinned 
the Swift XRT Photon Counting data with 250 counts 
per bin for the times up to 2 Ms after the burst and 100 
counts per bin for the times thereafter. Spectral analy- 
ses were performed for the times 300 - 800 ks and T > 
1 Ms after the burst. Source photons were collected in 
a circle with r = 23'.'5 and background photons with 
r = 96 with grade selection 0-12. The response matrix 
swxpc0tol2s0-20010101v010.rmf was used. The spectra 
were rebinned with 20 counts per bin for the 300-800 ks 
after the burst spectrum and 15 counts per bin for the 
spectrum with T > 1 Ms. The spectra were analyzed 
with XSPEC version 12.4.0x (jAinaud 1991). To search 
for changes in the X-ray spectrum we applied a hard- 
ness ratio studj0 segment by segment. Because of the 
low-number statistics in some of the later segments of 
the Swift XRT and all Chandra observations, we applied 
Bayesian sta tistics to determi ne the hardness ratios as 
described bv lPark et all ()2006D . 

3. RESULTS 
3.1. Temporal breaks 

The late-time X-ray light curve including the five 
Chandra pointings is shown in Figure[T] We ignore the 
first day of the Swift observation because it is not rel- 
evant for the study of the late-time light curve. The 
ea rly light curve and a detailed discussion can be found 
in IGrupe et al.l ([2007h . 

The late-time light curve (Figure[T]) was fitted by sev- 
eral power law and multiply-broken power law models 
as listed in Tabled Fitting the light curve with the de- 
cay sl ope aj3 fixed to 1.32, as reported bv IGrupe et al.l 
()2007[ ) from the Swift data, gives a very poor result {x^/v 
= 897/46). The fit can be improved by leaving the decay 
slope as a free parameter (Tabled model 2). This results 
in a single decay slope of as = 1.45±0.01, but the light 
curve still deviates significantly at later times from this 
slope, resulting in an unacceptable ^ jv = 400/45. A 
broken power law fit to the entire late-time light curve 
(model 3) reveals a break a t about 2 M s after the burst; 
in contrast to the result of IGrupe et al.l ()2007D . in which 
we could fit the late-time Swift data with just one decay 
slope, the addition of the 2007 Chandra data requires a 
break in the Swift data. The late decay slope, — 1.85, 
is driven by the last two Chandra observations, while the 
is also strongly affected by two very high data points 
at ~ 2 Ms and ~ 5 Ms. Making the assumption that 
these two high points are late-time X-ray flares unrelated 
to the afterglow of the external shock, we removed them 
from further fits (see models 4 and 5). We then fit the 
data between 100 ks and 30 Ms with a broken power law 
(model 6), obtaining a break time of Tbreak,3 = l-01^a22 

Ms and slopes of = l-32lo.o5 and = lMto,le- 
This fit is plotted as the dashed line in Figure[T] 

* The hardness ratio is defined as HR=(H-S)/(H-|-S), where S 
and H are the observed counts in the 0.3-2.0 and 2.0-10.0 keV 
energy bands, respectively. 

^ Wc fol lowed the n otation as defined in INousek et al.l 1)20061 ') 
and Zhang ct al.l II2006I ) where 03 is the steep decay slope after the 
plateau phase. 
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The last two Chandra observations deviate from this 
fit, suggesting a break at about a year after the burst. 
Because these last two points have very few counts (and 
consequently large uncertainties), a broken power law fit 
to the late-time light curve cannot constrain either the 
break time or the late-time decay slope as unless at least 
one parameter is fixed. We therefore approached the 
question of a final break in steps. A single power law 
fit to the light curve for T > 1.2 Ms (Table[2 model 7) 
gives a4 = 1.68±0.08 and x^/t^ = 12/15. Although this 
is already an acceptable fit, we investigated the possi- 
bility of a late-t hne break which is expected from GRB 
theory (c.f. e.g. IZhang et al.ll2006l: [Mfezaros II2006D . At 
first we fitted the light curve for 1.2 Ms < T < 35 Ms 
with a single power law (model 8) which results in 
= IMto^ and xVi^ = 6/13. We then fitted a broken 
power law model to the entire light curve with T > 1.2 
Ms with a4 fixed at 1.61 (the best- fit result when the last 
two Chandra observations are excluded; model 8 in Ta- 
blc[2]) to determine whether the data require a very late 
break in the light curve slope. This fit (TablcjU model 
9) gives rbrcak,4 = 41.3^51 Ms and 05 = 4.65^1% with 
X^/i^ = 6/14. Figure|3] displays the contour plot between 
the final break time and the final slope. It shows that 
they are still not well-constrained. Although the best-fit 
break time is 41 Ms (2007 November), a break as early 
as ^ 26 Ms, with a late-time decay slope of as = 2.5, is 
consistent with the data at the la level. 

In addition to the broken power law fits with a sharp 
break, the late-time light curve was also fitted by the 
smoothed double - broken power law model defined by 
IBeuermann et all (jl999D . Here we found decay slopes 
a3=1.20±0.07, a4=1.70±0.05, and as=2.41±0.26 with 
break times at 1.09±1.01 Ms and 20 Ms (fixed). The 
smooth parameter is fixed to 3.0 and 2.0 for the breaks 
at about 1 Ms and 20 Ms, respectively. This results in 
an acceptable fit with x^/z^ = 42/32. Possible interpre- 
tations of these temporal breaks are discussed in ^ 

3.2. Spectral variations 

Temporal bre aks are of t en as s ociated with spectral 
breaks (e.g ISari et all 119981: IMeszaros eraLl 119981: 
IZhang et al.l \20m . Figure^] displays the Swift XRT 
count rate and hardness ratio light curves for the interval 
between 100 ks and 5 Ms after the burst. The hardness 
ratios are plotted segment by segment. While the hard- 
ness ratios before the break at 1 Ms after the burst are 
of order HR ^ 0.3, after the break the spectrum hardens 
to HR ^ 0.45 with even harder values at later times. 

The spectrum before the 1 Ms break can be fitted with 
a single absorbed power law with N^i = (l-34to;25) x 10^^ 
cm~^ and an energy spectral slope /3x = 1.18±0.11 (j^ jv 
= 81/82). The spectrum after the 1 Ms break was also 
fitted by an absorbed single power law model. Leaving 
the absorption column density as a free parameter, how- 
ever, results in an increase of the column density, which 
does not seem plausible. Therefore wc fixed the absorp- 
tion column density to iVn = 1.34 x 10^"'^ cm~^, the value 
obtained before the break. This fit results in a slightly 
fiatter energy spectral slope /3x = 0.89±0.11. These val- 
ues were used in PIMMS to convert the Chandra ACIS-S 
count rates into the fluxes given in Tablc[T] and plotted 
in Figure[T] 



The Chandra data must be analyzed in the Poisson 
limit, complicating proper analysis of possible spectral 
variations at very l ate times. Using the Bayesian ap- 
proach described bv lPark et al.l (|2006D . we estimated the 
hardness ratios in the Chandra dateP^ and their uncer- 
tainties, both before and after the break at 38 Ms. We 
obtain mean values of HR=— 0.39 for the 2007 March- 
June data (before the final break; 38 counts total) and 
HR=— 0.80 for the very late data (after the final break; 8 
counts total), with 85% confidence limits of HR= —0.60 
to -0.17 and HR= -1.00 to -0.58, respectively. Al- 
though this is a suggestion of spectral softening across 
the final break, we cannot exclude (at the 85% confidence 
level) the possibility that the hardness ratio is constant. 
Note that due to the different energy bands and detector 
response matrices it is not possible to compare the Swift 
and Chandra hardness ratios directly. 

3.3. Comparison with other GRBs 

Even though GRB 060729 was one of the brightest 
bursts detected in X-rays, it is not the brightest one so far 
seen during the Swift mission. GRB 061121 was about 
2-3 tim es brighter when the Swift-XRT started observ- 
ing it (jPage et al.l 120071 ). but by a day after the burst 
it was an order of magnitude fainter than GRB 060729 
and was not detected after 2 Ms post-burst. The second 
brightest X-ray afterglow so far seen by Swift was GRB 
080319B (Racusin et al., 2008). Figure 11 displays the ob- 
served count rate light curves of GRBs 060729, 061121, 
and 080319B. Even though GRBs 061121 and 080319B 
appear to be much brighter than GRB 060729 until about 
20 ks after the trigger, the long plateau phase in GRB 
060729 makes it the brightest X-ray afterglow about half 
a day after the trigger. GRB 061121 already displays a 
break from the plateau to the second steep decay phase 
at 2.2 ks after the burst followed by an even steeper decay 
at 30 ks with a decay slope of a = 1.5 (Page et al., 2007). 
This earlier break compared to GRB 060729, which broke 
at about 60 ks after the burst, and the steeper decay 
slope of 1.5 compared to 1.3 in GRB 060729 made GRB 
061121 disappear much faster than GRB 060729. The 
"naked-eye" burst 080319B on the other hand docs not 
even show a noticeable plateau phase and decays rather 
quickly with late-time decay slopes of ax= 1-17 and 2.61 
before and after the late-time jet break at T = 9.5 x 10^ 
s. 

We can ask how the intrinsic rest-frame 2-10 keV lumi- 
nosity light curve of GRB 060729 compare s to the bursts 
shown in Figure 1 in iNousek et al.l (12006') . Figure[5] dis- 
plays GRBs 060729, 061121, an d 080319B with mos t of 
the bursts shown in Figure 1 in INousek et al.l ()2006[ ) in 
the rest-frame. The plot shows that during the plateau 
phase GRB 060729 was not that luminous. As a matter 
of fact it was a factor of about 10 to 100 less luminous 
than bursts such as GRBs 061121 or 080319B. Never- 
theless after a few days in the rest-frame, GRB 060729 
becomes the most luminous X-ray afterglow in the 2-10 
keV band. 

Compared with other GRBs shown in Figure[5l the to- 
tal energy output in the apparent rest-frame 2-10 keV 

1" We defined the Chandra hardness ratio as HR=(H-S)/(H+S), 
where S and H are the counts in the 0.5-2.0 and 2.0-8.0 keV bands, 
respectively. 
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band of £'2-iokoV = 7 x 10^^ ^rgs makes GRB 060729 
one of the most energetic X-ray afterglows ever detected. 
In the 2-10 keV band, only GRBs 061121 and 080319B 
appear to be more energetic with £'2-iokoV = 1 x lO^'' 
ergs and 2 x 10^^ ergs, respectively. However, if we 
attempt to correct for the jet opening angle, the pic- 
ture changes. With jet opening angles of 4° and 0.4° 
as inferred for GRBs 061121 and 08031 9B, respectively 
(|Page et all 120071: iRacusin et all |2008[ ). the beaming- 
corrected energies are 2.4 x 10^° and 4.0 x 10^" ergs, 
respectively. For GRB 060729, however, assmiiing a jet 
half opening angle of 14° (see next section) the beaming 
corrected energy is still 2.1 x 10^^ ergs in the rest frame 
2 - 10 keV band. 

4. DISCUSSION 

While our original Swift XRT light curve after the 
plateau ph ase seems to be con sistent with a single power 
law decay (jGrupe et al.ll2007[ ). the Chandra data make 
it apparent that there had to be a break at about 1 Ms 
after the burst. We showed above that this break coin- 
cides with a significant change in the X-ray spectral slope 
from I3x— 1.2 before an d 0.9 after that bre ak. Unfortu- 
nately, as mentioned in lGrupe et al.l (|2007f ). Swift could 
only follow the afterglow in the UVOT with the Wl filter 
up to a month after the burst due to bright stars in the 
UVOT field which caused some scatter in the UVOT and 
an enhanced background at the position of GRB 060729. 

According to Model 9 (excluding the two X-ray flare 
points at 2 Ms and 5 Ms) in Table 2, the late time X- 
ray afterglow is described by Fi.^^ cx t~°'^v^^\ where 
(1) = 1.32;^ °^, ^3 = 1.18 ± 0.11 for 105 ks < t < 
Tbrcak,3 = l-Oli^'l Ms; (2) a4 = 1.6lt^;i^, = 0.89 ± 

0.11 for Tbreak.3 < t < Tb,eak,4 - ^-l^t'^fl X 10^ S; (3) 
«5 = 4.65t2-0^ for t > rbrcak,4. 

The closure relation for 0.1 Ms < t < 1.0 Ms re- 
sults in as — 1.5/33 = —0.46 ± 0.17, consistent with 
both the ISM and wind models with theoretical expec- 
tation of as — 1.5/33 = —0.5 if vx » maxj^^c, '^m}, 
so as = (3p — 2)/4 and /S^ ^ p/2. The power law in- 
dex of the electron energy distribution is p ~ 2.43^'lQ Qy 
(2.36±0.22), derived from a3(/33). The closure relation 
for 1.0 Ms < t < 41 Ms results in 04 - 1.5^4 = 0.28;5j;i^, 
moderately consistent with the theoretical expectation of 
a4 — 1.5/^4 = 0.5 if the environment is a free win cQand 
the spectral regime is < < Vc, so = (3p — l)/4 
and /34 = (p— 1)/2. In other words, the break at ~ 1.0 Ms 
can be interpreted as a cooling break in the wind medium 
scenario in which the cooling frequency crosses the X- 
ray band. The power law index of the electron energy 
distribution derived from the value of an (the spectral 
index of this epoch has a relatively large uncertainty) is 
p = 2.48lo!jg, quite consistent with the value derived 
during the previous epoch. Therefore, a wind model is 
preferred from the observations before t = A\ Ms, break- 
ing the degeneracy of the wind and ISM models which 

For the ISM model, the closure relation is 04 — 1.5/34 = —0.5 
or 0.0, which can be excluded at the 5a- and 2a confidence level, 
respectively. Even if we adopt an — 1.5/34 = 0, the derived 
p ~ 3.15^q' JI/2.78 ± 0.22 from the temporal/spectral index is in- 
consistent with the p value derived with 03 and at the earlier 
stage. 



are b oth consistent with the earlier data (jGrupe et al.l 
I2OOI . 

The X-ray afterglow light curve for t > 2.75 x 10^ s 
excluding the two flares can be also well fitted by a sin- 
gle smoothed broken power law: (1) as = 1.31 ± 0.05, 
= 1.18 ± 0.11 for 275 ks < i < Tbrcak 3 = 2.43 ± 0.79 
Ms; (2) a4 = 1.96±0.09, /34 = 0.89±0.11 for Tb,.oak,3 < t 
(xt/doi = 42/32, smoothness parameter s = 3). The 
steepening of the decay at t = 2.4 Ms, if not due to the 
cooling frequency passing through the observing band 
(which results — = 0.25, inconsistent with the 
observation), should originate from the post-jet-break 
evolution. The change of the temporal decay index, 
a4 — Q3 ~ 0.65 ± 0.10, is consistent with a non-spreading 
jet break in a wind model with Aa = 0.5 (if the jet 
has significant sideways expansion, then the value of a^ 
should be equal to p ^ 2.4 — 2.5). However, the spectral 
hardening around the break time can not be well inter- 
preted in such a model. Furthermore, the transition of 
a jet break in the wind medium usually takes two orders 
of magnitude in time (Kumar & Panaitescu 2000), which 
is inconsistent with the observation of GRB 060729. In 
conclusion, the temporal break at i ~ 1 — 3 Ms is prob- 
ably not a jet break. 

There are two possible interpretations for the last ten- 
tative light curve break at t = 41 Ms, as follows: a jet 
break, or a spectral break in a spherical model. We dis- 
cuss these next, after which we consider the implications 
of the long plateau phase. 

4.1. Jet model 

The jet -|- wind model predicts a^ = p ^ 2.4 — 2.5 
for a sideways expanding jet or — a^ + 0.5 ~ 2.1 for 
a non-sideways expanding jet. The value of the model- 
predicted temporal index after the jet break thus can not 
be excluded at even the la confidence level (see Fig. [3]). 
Recently Zhang & MacFadyen (2009) have performed 
two dimensional simulations and calculations of GRB af- 
terglow hydrodynamics and emission. They showed that 
the sideways expansion of GRB jets can be neglected dur- 
ing the relativistic phase (for the sideways expansion of 
GRB jets see also Kumar & Granot 2003, Granot & Ku- 
mar 2003) and that the change in decay slope was larger 
than predicted analytically. Their results may further 
alleviate the above problem of the relatively shallow the- 
oretical slope compared with the observed steep slope. 
However, the spectral softening revealed by the hardness 
ratio evolution around this break time somewhat disfa- 
vors the jet break interpretation. 

If we interpret this break as a jet break, then a 
half-opening angle of the jet can be inferred (e.g., 
iChevalier fc Lill2000[) . Under the thin-shell approxima- 
tion for the post-shocked fluid of a relativistic blast wave, 
the conservation of energy (neglecting the initial baryon 
loading in the fireball) reads 

Ekdso = A/swF^c^ = constant, (1) 

where F is the bulk Lorentz factor of the downstream 
fluid just behind the shock front and F^ — a/2F is the 
Lorentz factor of the shock. The swept-up circum-burst 
mass by the blast wave is 

= / 'inr'^n{r)mpdr = -AR^-^nip, (2) 

Jo 6 — k 
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where the environmental density n(r) = Ar ^ . and A = 
no for the ISM case (fc = 0) and A = 3 x 10^^ cm^M* 
for the steUar wind case (fc — 2). 

The evolution of the shock radius R measured in the 
observer's frame is dR ~ 2r^crfi/(l + z), therefore, 



R 



{3 - k){4 - k)EkMot 



and 



(1 + z)T:AmpC 
(3 - k)Ek,iso 



l/(4-fe) 



1/2 



(3) 



(4) 



Inserting the inferred values of the physical parameters, 
we have 



K- L.ZX iU -C'fc4so,54"0,^l 



t 



41 Ms 



1/4 



cm. 



r-i.oi?,aso,54«o,-i Utm^ 



-3/8 



(5) 



(6) 



for the ISM case, and 



i? = 1.0 X 10^°i?^/L^,,A-L/,^ 



t 



41 Ms 



1/2 



cm, (7) 



r = 4.0i?^/t_54A; 



1/4 
,-1 



t 



41 Ms 



fjct 



= 14.0°£;, 



-1/4 4I/4 
fc,iso.54^*,-l 



t 



-1/4 



1/4 



41 Ms 



(8) 



(9) 



for the wind case. The values of Ek^iso and A^ adopted 
here can be found below. We adopt the convention of 
Qx = in cgs units. Unless Uq <^ 0.1 cm~'^, the 

jet has already decelerated to be non-relativistic (F ~ 2) 
while the fact that there was no jet break before t = A.lx 
10^ s argues against the jet in an ISM medium. In other 
words, the outflow of GRB 060729 is likely spherical if 
the circum-burst medium is ISM. However, in this way 
the temporal break at rbreak,4 can not be explained with 
the hydrodynamic/geometry effect. 

At such a late time, the jet may also enter the non- 
relativistic phase in a stellar wind medium and the hy- 
drodynamics is described by the self-similar Sedov - von 
Neumann - Ta ylor solution. T he non-relativistic transi- 
tion timfH is |Waxmanll2004[) 



isNT = 2.7Ejct,52A^ ^ yr. 



(10) 



where Ejot — Ek,iso(^fct/ ^ is the beaming-corrected 
kinetic energy of the outflow. The non-relativistic 
transition predicts a flattenin g in the light curve 
(|Huang fc Cheng|[200irZhang fc MacFadvenll2009l ). The 
theoretical temporal slope is as = {7p— 5)/6 ~ 2.0 — 2.1 
for p ~ 2.4 — 2.5, which is not inconsistent with the ob- 
servations if the large error bars of the observed slope are 
considered. However, the trend of steepening after the 

1^ Zhang & MacFadyen (2009) suggested that the beginning 
time of the SNT self-similar evolution is ~ 5tNR, where ^ 
200-Ei 54^47^ years is the time when a GRB jet with negligible 
sideways expansion becomes non-relativistic. 



break contradicts with the trend of flattening when the 
jet becomes non-relativistic. 

The jet is still inside the free wind bubble ~ 450 days 
after the burst in the observer's frame, indicating that 
the termination shock radius of the wind bubble is larger 
than ^ 32 pc. The size of a GRB progenitor star wind 
bubble depends on the density and pressure of exter- 
nal interstel lar m edium , and the mass l oss history of the 
progenitor (iGar cia-Seg ura et al.l Il996t iDai fc Wul 120031: 
iChevalier et al.ll2004ll . The large size of the wind bub- 
ble surrounding the progenitor of GRB 060729 suggests 
(1) a fast and tenuous wind during the RSG phase prior 
to the Wolf-Rayet phase, (2) the life time of the Wolf- 
Rayet stage is relatively long, and (3) the density and 
pressure of the external interstellar medium should be 
low. Therefore, from modeling the afterglow, the pro- 
genitor of GRB 060729 is not likely to be located in a 
giant molecular cloud or a star-burst environment. Fur- 
ther deep optical observation of the GRB site and its 
host galaxy may help to test this prediction. 

4.2. Spherical model 

We assume that the initial shock-accelerated electrons 
at such a late time have a broken power law distribution. 



dN, 



die I 7h 7e 



7e 



7m<7e<7b, 

-P2 ry, ^ ^ > 

, 76 ^ 7e, 



where pi = p 2.4 is the low energy power law in- 
dex and P2 is the high energy power law index, 7f, is 
the break Lorentz factor of electrons. This assumption 
is much more realistic than the single power law as- 
sumption, especially for the late time when the shock 
Lorentz factor decreases to the order of unity /the shock 
is no longer ultra- relativistic (e.g., Hededal et al. 2004; 
Niemiec & Ostrowski 2006; Spitkovsky 2008). For sim- 
plicity, we assume Ri, = ib/lm remains a constant in 
time (Li fc Chevalier 2001). In this scenario, tc [vc = vx) 
is equal to 1.3 x 10^ s, th {vb = vx, i^t is the typical syn- 
chrotron frequency of 7b electrons) is equal to 4.1 x 10'' 
s. The last steep decay segment can be described with 
"5 = (3?32 - l)/4 and ^5 = (p2 - l)/2. The high energy 
power index is therefore derived to be p2 = 6.53l^'yg, 
and the inferred spectral index is /Ss = 2.771q gp. From 
tm < 10^ s and t;, = 4.1 x 10*" s, we constrain the param- 
eter Rb > (4.1 X 10^/10^)3/4 ^ 91. The above pi ~ 2.4, 
P2 ^ 6.5 and Rb > 91 are quite similar to those derived 
for GRBs 991208 and 000301C in Li fc Chevaher (2001). 

The synchrotron emission from a spherical relativistic 
blast wave can be described by (e.g.. Chevalier fc Li 2000; 
Wu et al. 2005) 



iy„. - 3.7 X 1012^2 _^^^,V2_ 



p.1/2 -3/2 

2.5^fc,iso,53''5 



(12) 



= 3.4 X lO^e-ig.^e^i/l^S^/^, 3^^24/2 Hz, (13) 



^^^^B.~2.5^k,iso,53^*''5 



1.5^ 

,1/2 



— 1/2 

A^t^ niJy, 



(14) 



where eg, £3 are the energy equipartition fractions of 
electrons and magnetic field respectively, A^, is the stellar 
wind parameter. We have considered the cooling effect 
on Vc by synchrotron-self-Compton scattering processes. 
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The crossing time t,„ in the optical band {vopt ^ 10^^ 
Hz) by I'm no later than t = 10^ s gives 

4-0.5eB,-2.5-Bfc,iso,53 < 7.3 X 10^ (15) 

which is easily satisfied, while the crossing time tc in the 
X-ray band {hi>x ^ 2 keV) by Vc gives 

^e.-o.^e'JU^^kMo^sAl = 2.24 X IQ-^. (16) 

The flux density at frequency i^c, F^^, should be equal 
to the observed 2 keV X-ray flux density of F2koV ^ 
3.5 X 10~^ mJy at tc ~ 10^ s, which reads 

^:.-L^b%.M:LM:'-' - 8.75 X 10^. (17) 

The above afterglow model parameters can not be tightly 
constrained due to the lack of more conditions/equations. 
However, a large late-time isotropic kinetic energy of 
Ek.iso ^ 10^^ — 10^^ ergs and a small wind parameter 
of A, ~ 0.1 are favored if we assume Ce ^ 0.1 and 
eB ^ 3x 10~^ as inferred from fittings to other after- 
glows (e.g., Panaitescu & Kumar 2001; 2002). Within 

this parameter set, Rb = 5.3 x lO^'ejio.i^B^-t s-^sT 130*54 
is much higher than those derived from GRBs 991208 
and 000301C by Li k Chevalier (2001). 

4.3. Energy Injection in the shallow decay phase 

A large Ek;iso (> 10^^ ergs) at very late time (after the 
plateau phase ends) can be understood if the early long 
plateau phase detected both in X-ray and optical bands 
from a few hundred seconds to ^ 6 x 10^ s since the 
trigger is interpreted as the period of eii ergy injection, 
as already discussed in lGrupe et al.l (|2007f ). Possible sce- 
narios for this enhance d energy injection are re freshed 
shocks as suggested by iRees fc M eszarosi (|1998[ ) or the 
continuous energy injection by the strong magnetic field 
of a newbor n fast-rotat i ng m agnet ar or black hole as 
sugges ted by iDai &: Lul (|1998[ ) and I Zhang &: Meszarosi 
()200lD . In the latter model, the temporal index q of the 
injected luminosity (L oc t~'^ , E-mj (x t^^'') can be derived 
from the observed values of the temporal and spectral in- 
dices during the plateau ph ase (e.g.. Tab l e 2 of Zhang et 
al. 2006). From Table 5 of iGrupe et all (|2007l) . we can 
obtain g = 2(a2-/32 + l)/(l + /32) = -0.037± 0.101 with 
a2 = 0.14 ± 0.02 and = 1.18 ± 0.11 for the X-ray af- 
terglow plateau assuming i/m < Vc < vx- Such a q value 
is consistent with the theoretical expectation {q = 0) 
of the Poynting-flux dominated injection model. There- 
fore the isotropic kinetic energy of the shock increases 
linearly with time, Ek^iso cx; t. The prompt isotropic 
gamma-ray energy release in the 1 keV - 10 M eV band 
is E-y^iso = 1.6 X 10^2 grgs (jGrupe et all (|2007[ )'). Adopt- 
ing a typical GRB efficiency of 10% — 90%, the initial 
isotropic kinetic energy remaining in the afterglow shock 
is about the same order of Ej, 

iso; i.G. Ek i^Q i ~ 10 ergs. 
After the energy injection is finished, the shock energy 
Ek,iso is increased by a factor of ^ 100 compared to its 
initial value, so Ek^iso is the order of ^ 10^'^ ergs. The 
lack of a jet break up to t = 642 days results a lower 
limit of the half-opening angle of the GRB jet, i.e., 

0i > ^5°e;-LmA'JU (18) 



which corresponds to a beaming-corrected jet (double- 
sided) energy of 

£;jet > 3.4 X 10''eI/.1,^A:^J^ ergs. (19) 

This makes it one of the most energetic jets ever seen 
and is the reason why we were able to detect the X-ray 
afterglow still 642 days after the burst. Note that the 
energy of a fastest-rotating magnetar is Em ~ 5-^^^ ^ 
2 X lO^^XisPj^g^ ergs. A massive and fast-rotating black 
hole as the central engine is thus more likely to provide 
the required energy. 

5. CONCLUSION 

What makes the X-ray afterglow of GRB 060729 so re- 
markable is the fact that it was still detected even almost 
two years after the burst. This exceptional late-time de- 
tectability is related to three things: a) with an initial 
0.3 - 10.0 keV flux of almost 10~^ ergs cm~^ it was 
one of the brightest af terglows ever detected by Swift 

b) its flat decay phase (jNousek et al.ll200"6l : iZhang et al.l 
I2OO6D extended out to about 60 ks after the burst, and 

c) the decay slope after that break is about a = 1.3. 
Despite breaks at T ~ 1 Ms and T ^ 1 year the af- 
terglow was still detected by Chandra nearly two years 
after the burst. Bursts like GRBs 060614, 061121, or 
even 080319B (iMangano et al.l [20071: iPage et al. I [20071: 
IRacusin et al.l [20081 respectivelv) were even brighter in 
X-rays at about 100 s after the burst than GRB 060729, 
but their plateau phases are significantly shorter than 
that of GRB 060729. They therefore faded more rapidly 
than GRB 060729 at late times. 

Analysis and modeling of the X-ray afterglow of GRB 
060729 show that this burst happened in a tenuous wind. 
This is consistent with the coUapsar picture of long 
GRBs. During the early plateau phase, the energy in 
the external shock increased by two orders of magnitude. 
A reanalysis of the Swift XRT light curve together with 
the three detections by Chandra in 2007 reveal a tempo- 
ral break at ~ 1.3 Ms after the burst. The decay slope 
steepened from a ~ 1.32 to a = 1.61 around this break 
and the X-ray spectrum hardened in the meanwhile, in- 
dicating this break is a cooling break (the cooling fre- 
quency of synchrotron radiation crosses the X-ray band). 
There is another light curve break at ~ 1.3 year after the 
burst tentatively indicated by the last two Chandra de- 
tections. This break coincides with a possible spectral 
softening, suggesting that the break may be of spectral 
origin, though a hydrodynamic origin (jet break) is also 
possible. If due to a jet break, then the implied half- 
opening angle is dj ^ 14°. If due to a spectral break, 
such a spectral softening could be the result of a very 
steep power-law distribution of shock-accelerated elec- 
trons responsible for the synchrotron radiation. In this 
case, with no evidence for a jet break up to 642 days after 
the burst by Chandra the jet half-opening angle must be 
9j > 15° and the jet energy Ej > 3 x 10^^ erg. Such a 
large jet energy implies that the central engine must be 
a fast-rotating massive black hole, not a magnetar. 

Our Chandra observations presented here have shown 
again how important Chandra is for late-time obser- 
vations of GRB X-ray afterglows. Chandra has al- 
ready been essential for the detection and non-detection 
of jet breaks in the X-ray afterglows of the short- 
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duration GRBs 0507 24 and 05 1221 A (jGrupe et al.|[2006l: 
iBurrows et al.|[2006l respectively). 
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Time since Burst 



Fig. 1.— Swift-XRT and Chandra ACIS-S light curve of the X-ray afterglow of GRB 060729. The black crosses di splay the Sw i ft XRT 
data and the red triangles the Chandra data. The solid line displays the initial decay slope 03 =1.32 as reported bv lGrupe et al.l l|2007l ). 
the dashed line the decay slope post-break at 1.0 Ms after the burst with 04=1.61 (from model 6 in Table 2), and the dotted line the steep 
decay after the jet break at 41 Ms after the burst with as = 4.65. 
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Fig. 3. — Contour plot between the late-time break and deeay slope 05 of the fit to the late time light curve of GRB 060729 with T >1.2 
Ms after the burst, excluding the flares at 2 and 5 Ms (Model #9 in Table[2]l. The lines mark the 1, 2, and 3 o" confidence levels. 
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Fig. 4.— Comparison of the observed X-ray light curves of GRBs 060729, 061121, and GRB 080319B 
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Fig. 5. — Rest frame isotropic equivalent luminosity light curves of several Swift and pre-Swift bursts taken from [Nousek et al.l 1)2009 ) 
including GRBs 060729, 061121, and 080319B. 
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TABLE 1 

Chandra Observation log of GRB 060729 



ObsID 


T-start^ 


T-stopi 


T 2 

cxp 




.Pb.S-lO.OkcV* 


7567 
8541 


2007-03-16 11:39 
2007-03-17 13-02 


2007-03-16 19:57 
2007-03-17 14-54 


27690 
4701 


3 

•J-fJ-O.gi 


4 26+^-^^ 


7568 


2007-05-16 09:30 


2007-05-16 20:56 


40268 


-0.66 


2 90+^-"^ 


7569 


2007-06-30 06:19 


2007-06-30 23:48 


60400 


2 07+"-''^ 


9 9t;+U.V8 
^■^'-'-0.49 


9086 
9801 
9802 
9803 


2007-12-26 01:02 
2007-12-28 19:23 

2007- 12-29 18:57 

2008- 01-05 01:04 


2007-12-26 09:04 
2007-12-28 23:59 

2007- 12-30 01:13 

2008- 01-05 04:23 


27395 
15068 
20239 
9986 




u.Da_Q 2g 


0.75l«:i« 


9087 
9811 
9812 
9813 


2008-04-30 23:09 
2008-05-01 22:45 
2008-05-03 15:51 
2008-05-04 15:07 


2008-05-01 08:34 
2008-05-02 07:50 
2008-05-03 23:35 
2008-05-04 23:22 


32330 
30775 
26382 
27810 




0.28t«:i| 



^ Start and End times arc given in UT 

2 Observing time given in s 

3 Count rate in units of 10~* ACIS-S counts s'^ in tfie 0.5 - 8.0 koV band 

* Unabsorbcd 0.3 - 10.0 keV fiux in units of 10^^" ergs s^^ cm^^. The count rates were converted assuming A'^^f = 1.34 X 10^^ cm"^ and 
/3x = 0.89. 



TABLE 2 

Fits to the late-time {T> lO^s) Swift XRT and Chandra light curve of GRB 060729 



Model 


"3 


^tarcak,3 




-^tarcak,4 


"5 




1) powl fixed^ 

2) powl free^ 

3) bknpowl 

4) powl* 

5) bknpowl* 

6) bknpowl T < 30 Ms* 

7) powl T > 1.2 Ms* 

8) powl 1.2 Ms < T < 35 Ms* 

9) bknpowl T > 1.2 Ms*'^ 


1.32 (fixed) 

1.45±0.01 

1 ,9+0.03 
^■■^^-0.01 
1.46±0.04 
1 09+0.02 
^^■•-'^-0.05 

1 09+0.02 

^•■-'^-0.05 


1 90+0.25 
^•'^■'-0.40 

^^•'-'^^-0.22 


-, otr+O.lO 

-, 7Q+0.06 
^- '^-0.07 

J^-OJ^-0.06 

1.68±0.08 

^■"^-0.13 
1.61 (fixed) 


4i.3t*:? 


*-°^-1.34 


897/46 

400/45 

168/43 

160/41 

65/39 

60/37 

12/15 

6/13 

6/14 



^ Break time Tbreak a^re given in units of Ms 

2 Decay slope fixed to the value given in IGrupe et al.l {2003) 

Decay slope parameter left free to vary 
* Excluding the two flares at 2 Ms and 5 Ms 

^ Error bars of 05 and T^reak 2 a^re determined by keeping 04 fixed at the best-fit value for model 4 and assuming that ^ • follows 
the Gaussian distribution (see Fig. [Sjl. 



